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^ ' nance of traversable worniholes. In this paper, we explore the possibility that 

if such fields exist, they might be used to produce large negative energy fluxes 
and macroscopic violations of the generalized second law (GSL) of thermody- 

/\^ • namics. We find that it appears to be very easy to produce large magnitude 

negative energy fiuxes in flat spacetime. However we also find, somewhat sur- 
prisingly, that these same types of fluxes injected into a black hole do not 
produce violations of the GSL. This is true even in cases where the flux results 
in a decrease in the area of the horizon. We demonstrate that two effects are 
responsible for the rescue of the GSL: the acausal behavior of the horizon and 
the modiflcation of the usual black hole entropy formula by an additional term 
which depends on the scalar field. 
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1 Introduction 

The weak energy condition (WEC) states that [jl|: T^yU'^U" > 0, for all timelike vec- 
tors, where T^,y is the stress-energy tensor of matter, and U^ is an arbitrary timelike 
vector. By continuity, the condition holds for all null vectors as well. Physically the 
WEC implies that the energy density seen by all observers is non-negative. All exper- 
imentally observed forms of classical matter satisfy this condition. However, quantum 
matter fields can violate this and all the other standard energy conditions of general 
relativity 0. In at least some circumstances, quantum field theory imposes certain 
restrictions on the degree of energy condition breakdown, in the form of what have 
come to be known as "quantum inequalities" . These typically restrict the magnitude 
and duration of negative energy densities and fluxes, and seem to limit the production 
of gross macroscopic effects. The quantum inequalities put severe restrictions on the 
realizability of traversable wormholes and warp drives 0, ^, ^, ||, |^ . 

It is often assumed that all classical fields will obey the WEC. This is, however, 
not always the case. The energy density for a classical non-minimally coupled scalar 
field can in fact be negative. To our knowledge, the associated stress tensor was 
first written down for the case of conformal coupling by Chernikov and Tagirov [Q. 
The advantages of conformal coupling were emphasized by Callen, Coleman, and 



Jackiw [0 for flat spacetime quantum field theory, and further discussed by Parker |10| 



in curved spacetime. Bekenstein [[T^, noted that the conformal scalar field can 
violate the WEC and other energy conditions, and took advantage of this property 
to construct non-singular cosmological models. The violation of energy conditions 



by non-minimal scalars was further discussed by many authors, including Deser [O 



Flanagan and Wald WM , and in more detail recently by Barcelo and Visser |TB|, ITBI, ITT 



These features have enabled Barcelo and Visser to construct macroscopic wormhole 
solutions in which non-minimally coupled scalar fields violate the averaged null energy 
condition and serve as the "exotic matter" source required for wormhole maintenance. 
Because such fields are classical, they are not subject to the quantum inequality 
restrictions on the magnitude and duration of negative energy. 

This situation is rather unsettling, since the door is then opened for all sorts 
of bizarre effects. One of the most disturbing of these is a potential violation of the 
second law of thermodynamics by creating fluxes of negative energy. For example, one 
might shine such a flux at a hot object and decrease its entropy. If the radiation field 
has zero entropy, then the second law will be violated. In the case of negative energy 
fluxes produced by quantum fields, the quantum inequahties appear to prevent such 
large scale breakdowns of the second law [jl8|, |19[. Indeed in some cases, such as the 



Hawking evaporation of black holes, negative energy is required for the consistency 
of the unification of the laws of black hole physics and the laws of thermodynamics. 
The generalized second law (GSL) states |^ that the sum of the entropy of a black 



hole and that of any surrounding matter can never decrease, so 

^Stotal = ^S + ASmatter > (1) 

where here S is taken to be the entropy of the black hole, which is proportional to 



the area of its event horizon. In the Hawking evaporation process, although the area 
and thus the entropy of the black hole decrease, this is more than compensated for 
by the entropy of the emitted thermal Hawking radiation. 

In this paper we demonstrate that large, transient negative energy fluxes can be 
produced quite easily with classical free massless non-minimally coupled scalar fields, 
even in flat spacetime. Such fluxes appear to have magnitudes large enough to violate 
the second law for arbitrary lengths of time. However, due to possible uncertainties 
of how such fields might interact with ordinary matter, one could conceivably argue 
that perhaps the second law might still hold. In an attempt to circumvent the latter 
possibility, we examine a classical non-minimally coupled scalar field on a black hole 
background. We show that although the integrated energy flux injected into the 
hole is positive, it can be made temporarily negative, and that these periods can be 
made arbitrarily long. This raises the possibility that one might be able to achieve 
temporary but large, and in principle measurable, violations of the GSL. We give a 
proof that the GSL is in fact always satisfied, even though the horizon area of the 
black hole can temporarily decrease. A physical understanding of how this happens 
involves the subtle interplay of two effects - the acausal nature of the horizon and 
an additional term in the black hole entropy which depends on the scalar field. Two 
illustrative numerical examples are presented which demonstrate that a consideration 
of both effects is crucial to the preservation of the GSL. We conclude with a summary 
of our results and a discussion of open questions. 

In this paper, we use the MTW metric signature and sign conventions, and work 
in units where G = h = c = 1. 



2 Negative Energy Fluxes 

The Einstein equations for a generically coupled scalar field can be obtained by 
varying the Einstein-Hilbert action: 

S=^j d'^V^ R + J d^x^J^ [-]^g^^^ d^<l> 9,0 - \/(0) -\iR 02^ . (2) 
The resulting stress-energy tensor, T^j,, for the scalar field has the form 

T^u = V^0 V,0 - -(7^.(V0)2 - g^, \/(0) 

H [G^u 0' - 2 V^(0 V.0) + 2 g^, V\0 Va0)] , (3) 

where G^jy is the Einstein tensor, and V{(f)) is the scalar potential. If all the depen- 
dence on G^jy is grouped on the left hand side of the Einstein equations, then we can 
rewrite them by using an effective energy-momentum tensor given by 



T, 



"" 1 - 877^02 

-2e[v^(0V,0)- g^,V 



V^0 V.0 - -(?^.(V0)2 - g^u \/(0) 



(4) 



This is the expression used by Barcelo and Visser for their analysis of energy condi- 
tion violations. The Einstein equations now read G^y = SttT^^. Thus a constraint 
on this effective stress-energy tensor is translated directly into a constraint on the 
spacetime curvature. Note that T^,y will become singular as (1 — Svr^^^) — >• 0, unless 
the numerator in Eq. (^) vanishes at least as rapidly. In this paper, we will make 
the restriction that (1 — Svr^^^) > everywhere. The wormhole solutions of Barcelo 
and Visser have the feature that (1 — 87r^0^) can be negative, which implies either 
super-Planckian values for or extremely large values of the coupling parameter ^. 
However, we will show that even with the restriction that (1 — 87r^0^) > 0, one can 
produce negative energy fluxes of large magnitude. 

2.1 Fluxes in Flat Spacetime 

In the remainder of this section we will set V{(j)) = 0. Let us first demonstrate that 
the effective stress-tensor, Eq. (1), can lead to potential problems with the second 
law even in flat spacetime. We choose the simple case of waves travelling only in the 
positive a;-direction, where 

= 0(t - x) . (5) 

In this case, the energy density T**, and the flux T*^' = F^., are equal, and so the 
energy density will also be negative when the flux is negative. (To distinguish the 
case of a positive flux in the minus x-direction from a true negative energy flux, we 
count the flux as negative only when the energy density is simultaneously negative.) 
We will assume in this subsection that 

87re0^ < 1 . (6) 

With the above choices, the energy flux in fiat spacetime reduces to 

F, = r- = -T,. = [(1 - 20(0')' - 2^001 = (00' - a<py , (7) 

where we have used the fact that if has the form of Eq. (^ , 

dt(f) = -(9^.0 = 0' 



dxdi 



l" 



X'^t'i 



Let to be a point at which 0' has a nonzero extremum, so that 0'(to) = 0, but 7^ 
and 0" 7^ 0. Then, we can make F^; < in a neighborhood of to by arranging to have 
-2^00" < at to. 

As an example, choose ^ > and 

(f){t- x) = Bsmiu{t- x) +b, (9) 

where u = k = k^ > 0, and we assume that i? > 0, 6 > 0. Require that 

^»(5 + &)', (10) 



which will satisfy our condition Eq. (^). The instantaneous flux, when cosuj{t—x) = 
and sinci;(t — x) = —1, is 

F. = -2^5Vf4-lV (11) 



Let B^ = (3/i, with /3 < 1 and 6/^ - 1 ~ 0(1), which will satisfy Eq. (^. Then 

F.^ -2/3^2^-/3^, (12) 

where T = 2tt/uj is the period of oscillation. This expression is in Planck units. We 
can convert to conventional units by recalling that the Planck unit of flux is 

i^P = ^- 1.6x10^°^^, (13) 

l^ tp cm^ s 

where rup ^ 2.2 x 10~^g, Ip f» 1.6 x 10~^^cm, and tp ^ 5.4 x lO^'^^s are the Planck 
mass, Planck length and Planck time, respectively. Then Eq. ( |T^ ) becomes 

F. ^ -Sn'PFpS^^ -3.7 x 10^^ ^ /? (]^)\ (14) 

i ^ cm"' s V i / 

If, for example, we set T = Is, 6 = 2i? and f3 = 0.01, we have a negative energy flux 
of magnitude 

|F,|^3.7xl0i^-^, (15) 

cm"' s 

which has a duration of the order of a few tenths of a second. By ordinary standards, 

this is an enormous amount of negative energy, which has the potential to cause 

dramatic effects. 

Note that Eq. (^ shows that the flux can be expressed as a positive quantity plus 
a total derivative. This means that if and its first derivative vanish in both the 
past and the future, then the time-integrated flux is positive. As the above example 
illustrates, however, the transient flux can be made negative for an arbitrarily long 
time. 

In order to argue that this flux can be used to achieve a violation of the second 
law, we need to specify how the scalar field interacts with ordinary matter. It might 
be argued that if this interaction is sufficiently weak, then second law violations might 
be avoided. Without a detailed model of a system to detect the negative flux, this 
question may be hard to answer definitively. Our purpose here was to demonstrate 
that one can already achieve negative energy fluxes of large magnitude, which have 
the possibility of violating the second law, even in flat spacetime. In the next section 
we turn our attention to a relatively unambiguous energy flux detector - a black hole. 

2.2 Fluxes in Curved Spacetime 

In this section, we will consider the absorption of a classical scalar negative energy 
flux by a black hole. In Schwarzschild spacetime, from Eq. (^, a flux in the radial 



direction is given by 



^<0' 



:i - 200 



^t(P,r ■ 



-2e 



C. 



>,tr 



2C 



(16) 



where C = 1 — 2M/r, and C,r = IMjr^ . If we redefine the radial coordinate in the 
standard way, using r* = J C~^dr, we can express the right hand side as 



t; 



1 



1 - 2O0,t0,r* -2q 00,tr* -2'^"' ^^'^ 



1 - 8<02 

Consider an ingoing wave of the form 

where v = t + r* and [6, (f) are angular coordinates. Then 



(17) 



and 



T/ 



^7re02 



(1 -20(0,.)'- 2^(00,.. -^C„ 



(19) 



(20) 



The mass change of the black hole due to the absorption of the flux can be 
calculated from 

E = JT.^lS^dT.^ (21) 

where (3^^ is the timelike Killing vector and dH^ is the area element of the three- 
dimensional hypersurface defined by the black hole's horizon. Here we are assuming 
that the fractional change in the black hole's mass, M , is small over the time scale 
M, the light travel time across the hole. Therefore the metric is approximately 
Schwarzschild and hence has a timelike Killing vector. By energy conservation, the 
time rate of change of the mass is given by [^ 



M 



dM 
"dt 



T/ r^dn , 



(22) 



On the horizon, T/' = T„„, so this becomes 



h -Lt 



M 



m' / dn {Tnr=2M 



4M^ 



dQ 



(- 



1 



l(l-8<02) 
+ a positive quantity , 



dQ 



ln(l-87r^0^ 



(1-20(0,.) - 
1 



AM 



In 1 - 8n^4>- 



'"" ' 2M 

2 



=2M 



(23) 



where we have used the fact that, on the horizon, C,r= 1/(2M). As an example, if 
we initially assume that ^ > and > 0, then when 0,„ = and 0,^,^, > 0, we can 
have M < 0. 

Strictly speaking, r* -^ — oo at r = 2M . However, we can get around this problem 
by taking the effective boundary of the hole to be at r = 2M + e. That is, we draw a 
surface very close to the horizon, at a finite value of r*, and assume that any waves 
that pass through this surface also go into the black hole. Then we have 

■ dM dM 

dt dv ^ ^ 

As is the case for the flat spacetime flux, M can be written as the sum of a total 
derivative with respect to v and a positive quantity. Thus if and 0,^ vanish as 
V — > ±cxD, the net change in M is positive. However, one might worry that this still 
allows the possibility of macroscopic violations of the second law over long periods of 
time. 



3 Black Hole Entropy 

It has been shown by a number of authors |]22|, ^ ^ that in the presence of a classical 

field on the horizon, the entropy of a black hole is not simply 5* = A/4, where A is the 

area of the event horizon, but is modified by the presence of additional terms. Visser 

23[| has referred to these and similar objects as "dirty" black holes. For a classical 



non-minimally coupled scalar field, the Lagrangian density is given by 

^ = -^(^'^'V^0V.0 + ei?0' + W))- (25) 

In the case where the Lagrangian is an arbitrary function of the scalar Ricci curvature 
R, the entropy of a stationary black hole is given by (see Eq. (59) of Ref. p3| ) 



where C is the (Euclidean) Lagrangian density for the scalar field, and where we 
have set h = k =Boltzmann's constant= 1. The integration runs over a spacelike 
cross section of the horizon H. In this section, our arguments are not restricted to 
Schwarzschild black holes. 

In our case, dC/dR = —(1/2)^0^, so 



TrsV29d X = -- 
H oR 2 



2^rf\ = -^^(0^)A, (27) 



where 



^') = ^ U'V^d'x (28) 



A 



is the average of over the horizon. The resulting expression for the entropy of a 
stationary black hole is then simply 

^ = J (1 - Svre (0^)) . (29) 

One can interpret Eq. (B^) as 

A 



'^TGli,' '^°' 



where 



Ge// = (l-87re(0^))-^ (31) 

The scalar field effectively modifies the local value of Newton's constant. When ^ = 0, 
the case of minimal coupling, our expression reduces to the usual Bekenstein-Hawking 
entropy. The second term only contributes to the black hole's entropy when 7^ 
on the horizon. Here we are discussing the effect of the presence of the scalar field on 
the entropy of the black hole. Since the radiation field itself is a particular solution 
of a classical wave equation, its entropy Smatter is zero. 

Initially we assume that = on the horizon, so the initial entropy of the black 
hole is 5*0 = A/4. We then shine in a classical scalar field flux from infinity onto the 
horizon. During the time that 7^ on if , we get a contribution to the black hole 
entropy from the second term in Eq. (p9|). In order to count as a violation of the 
GSL, any decrease in entropy should be sustainable for a least a time of order M, 
in order to make a measurement of the horizon area. At a later time, the entropy 
might increase again due to oscillations in the sign of the flux, but we would still have 
achieved a measurable violation of the second law. 



3.1 Einstein vs Jordan Frames 

Here we address the issue of the existence of various conformal frames for scalar field 
theories. The form of the action of Eq. (0) is often called the Jordan frame. It is 
possible to perform a conformal metric transformation and field redefinition which 
convert the action into a form in which the ^-R0^ term is no longer present. To our 
knowledge, this was first done for the transformation relating a massless conformally 
coupled and a minimally coupled scalar field by Bekenstein |1T| , and later generalized 



to the case of arbitrary values of the coupling parameter by others. (See the excellent 
review article by Faraoni, Gunzig, and Nardone |2^ for details and further references.) 



Deser |T3| pointed out that this transformation can remove violations of the WEC 
The form of the theory in which ^ is transformed to zero is called the Einstein frame. 
The use of the word frame to describe these forms of the action is perhaps misleading; 
the transformation between the two is not a coordinate transformation, but rather a 
field redefinition which mixes the nature of the scalar and gravitational fields. 

This field redefinition technique has been extensively used by Jacobson, Kang, 
and Myers p6| to study black hole entropy, particularly in the context of higher 



curvature gravity theories. One could apply their technique to the massless non- 
minimally coupled scalar field case as well to argue that if the GSL is preserved 
in the Einstein frame, then it should also be preserved in the Jordan frame. The 
argument goes as follows ||2^. Each solution in the Jordan frame can be mapped, via 
the field redefinition technique, into a solution in the Einstein frame. In particular, 
the black hole entropy for a stationary state, as formulated in the Jordan frame, 
can be transformed to the entropy in the Einstein frame. The latter is just the 
usual Bekenstein-Hawking entropy, S = A/4, which we know is non-decreasing in the 
Einstein frame, provided that any matter absorbed by the hole satisfies the WEC 
Additionally, a non-minimally coupled scalar field (which violates the WEC) in the 
Jordan frame is transformed to a minimally coupled scalar field (which obeys the 
WEC) in the Einstein frame. Now suppose that in the Jordan frame the evolution of 
the black hole proceeds from one stationary state to another via a series of equilibrium 
states. Since each solution representing one of those states can be transformed into 
a solution in the Einstein frame, and since in the Einstein frame the entropy is a 
non-decreasing function, it should also be a non-decreasing function in the Jordan 
frame. Hence, the GSL should be satisfied in both frames. 

The question which now arises is whether the Jordan frame and the Einstein 
frame are physically equivalent to one another. This issue has been widely debated 
in the literature and is reviewed in Ref. |2^. There is certainly a formal equivalence; 
a solution of the Einstein-Klein Gordon equations in one frame can be mapped to 
a solution in the other frame. Thus in the Jordan frame, one must have a function 
which is non- decreasing in time obtained by mapping the Einstein frame area to 
the Jordan frame. What is less clear is the physical interpretation that should be 
attached to this function, or what observers in the Jordan frame would actually see if 
negative energy is absorbed by a black hole. Furthermore, one is often interested in 
situations where the full theory is not given by Eq. (^ alone, but there are additional 
parts in the action describing other fields. The transformation to the Einstein frame 
would then require transforming these parts as well, potentially leading to a very 
cumbersome formulation of the theory. 

For these reasons, it is desirable to have a proof of the GSL formulated entirely 
in the Jordan frame. Such proofs have been given by Jacobson, Kang, and Myers 
p6[| for gravity theories with actions that are polynomials in the Ricci scalar, and 



by Kang [^ for the Brans-Dicke theory. Their methods would presumably also be 
applicable to the case of non-minimal scalar fields. In the next subsection, we will 
give a somewhat different proof of the GSL for this case. 

3.2 Proof of the Generalized Second Law 

In this section, we will show that the entropy defined by Eq. (p9|) is a non-decreasing 
function, even when negative energy is being injected into the black hole. Strictly 
speaking, this is the entropy of a stationary black hole. Let us now assume that the 
expression given by Eq. (pOl) is the black hole entropy even at each moment during 



the dynamical evolution of the horizon and examine its behavior. We assume that 
0^0 outside of a finite interval in A, where A is an afiine parameter on the horizon. 
In general, can vary over the horizon, so it is convenient to divide the horizon into 
n patches, where each patch is sufficiently small that there is no spatial variation in 
through the interval in A during which 7^ 0. Thus S can be written as a discrete 
sum 



s = Y.^s, 



i=l 



where 



ASi = -AAi(l-87r^ 



(32) 



(33) 



Here AA, is the area of the z-th patch and 0j is the value of on that patch. Now 
we consider a particular patch, and drop the i subscript. The logarithmic derivative 
of AA with respect to A is the expansion, 6 (see, for example Ref. ||2^, p. 137): 

1 dAA rfln(AA) 



Thus we can write 



d_ 
dX 



AA dX 


dX ^''' 


ln(A5) =9- 


Stt^ d 
1 - Stt^ 02 dX 



(34) 



(35) 



We will take another derivative of this equation with respect to A and use the Ray- 
chaudhuri equation for the null generators of the horizon: 



de 

dX 



2 
2 



2 2 



ijj 



o 



K^iyk k 



8n T^^k^'kr 



(36) 



where k'^ is the tangent vector to a generator, a is the shear and uj the vorticity. 
From Eq. (HI), we can write 



(37) 



T,.k^k^ ^_g^^^,pr[v,0v. 


^-2eV^(0V. 


These relations can be combined to yield 


'^An(AS)= ^e' c^2 ^2 8^ 

dX^ ^ ' 2 1-871^02 


r 327reV' 1 
l-87re02_ 



dX 



(38) 



In this calculation, we have used the facts that d/dX = k'^Vfj, and that, because /c^ is 
the tangent to a geodesic, k'^V ^k^ = 0. So long as (1 — 87r^02) > 0, all of the terms 
on the right hand side of Eq. (^) are less than or equal to zero. Hence 



d'' 



dX 



2 ln(A5) < , 



(39) 



10 



Now we wish to show that AS* is a monotonically increasing function of A. First 
note that In (AS*) is a monotonic function of AS", so either one being monotonically 
increasing implies that the other is also monotonically increasing. We assume that 
the scalar field vanishes as A — *■ cxo, so that AS* —>■ \AA. Furthermore, at late times 
AA approaches a constant value of Af/n, where Af is the final horizon area. Thus 

-^ In(A^) ^0 as A ^ oo . (40) 

uA 

If AS* were ever to be a decreasing function of A, we would have to have d[\n{AS)]/dX < 
at some finite value of A. However, this is impossible because the only way a twice 
different iable function can go from having a negative value of its derivative to a zero 
value is to pass through a region where the second derivative is positive, which is 
forbidden by Eq. (^). Thus we conclude that 

AaS>0 (41) 

everywhere. Finally, if the entropy contribution of each patch is nondecreasing, it 
follows that the entire entropy is also a nondecreasing function: 

This completes the proof of the generalized second law. Note that V{(j)) does not 
appear in Eq. P7|), and hence our proof is independent of any scalar self-coupling. 
The proof is not restricted to Schwarzschild black holes, and so would also hold for 
other, e. g., Kerr, black holes as well. Here we also did not need to assume that the 
changes in the mass of the black hole must be necessarily small. 

We point out that although S* is a non-decreasing function, the same is not nec- 
essarily true of the horizon area. If we take the second derivative of ln(Ay4), and use 
the Raychaudhuri equation we can write 



dHn{AA) l/d\n{AA)^'^ 



dX^ 2 V dX 

If T^uk^'k" > 0, then 



uj^-a'-SnT^^k^'k''. (43) 



d^HAA) 

-^^<0. (44) 

However, if the WEC is violated, so that T^^k^k" < 0, then the sign of d^ \n{AA)/dX^ 
can be positive, and therefore the argument we made for the non-decrease of 5" does 
not work for A. We assumed that the scalar field vanishes outside of a finite interval 
in A, so we know that the entropy in the distant past is 5*0 = Aq/4 and in the distant 
future is Sf = Af/4, where Aq and Af are the initial and final total horizon areas, 
respectively. Since we proved that S" is a non-decreasing function, this implies that 
Af > Aq. The case Af = Aq corresponds to the trivial case that = everywhere, as 

11 



may be seen from the following argument. If Sf = Sq, the entropy must be constant. 
However, this can only occur if AS is constant for each patch, which in turn can only 
occur if = for all A. If AS* is constant, all of the terms on the right hand side of 
Eq. ( p8| ) must vanish, but the last of these terms can only do so for the trivial, (f) = 0, 
case. Of course, there is still the possibility of temporary decreases in A during the 
evolution. Indeed in the next section, we will exhibit specific examples where A can 
either increase or temporarily decrease, while S is always non- decreasing. Note that 
the fact that Af > Aq is a generalization of our result in the previous section. There 
we showed that the net change in the mass of a Schwarzschild black hole is positive 
in the case of small fractional changes. 

4 Numerical Examples 

In this section, we will construct some numerical examples to illustrate the effects of 
negative energy on a Schwarzschild black hole. We will restrict our examples to the 
case of spherically symmetric pulses which change the black hole's mass by a small 
fraction. 

4.1 General Procedure to Solve for 9 and A 

The relation between the expansion of the null generators, 6, and the area of the 

horizon. A, is given by 

rIA 

- = A<,(A), (45) 

Strictly speaking, in the general case, A is the area of an infinitesimal cross- sectional 
area element of a bundle of null geodesic generators. The expansion 6 measures 
the local rate of change of cross-sectional area as one moves along the generators. 
However, in this section, we assume spherically symmetric black holes and spherically 
symmetric energy distributions. Therefore the infinitesimal change in horizon area 
will be the same along each bundle of null generators, so in our case A can be taken 
to be the total horizon area. 

We also need the relationship of A to Schwarzschild-like coordinates. For a sta- 
tionary black hole, we can take (see, for example Ref. |^, p. 122) 



where k = 1/(4M) is the surface gravity of the black hole and M is its mass. Here 

V = e""" (47) 

is the Kruskal advanced time. Thus we have that 

A = ft- V = ft-i e"^ (48) 

12 



We will assume that T^^(f) represents ingoing radiation which changes the black 
hole's mass by only a small fractional amount, |AM| <C M. We can then take k 
to be constant to lowest order. If we change the independent variable from A to 
V = t + r* , then 

and 

1 p^pK fc ^ e J-vv ■ \y^) 

For spherically symmetric pulses, the shear and vorticity vanish, and the Ray- 
chaudhuri equation, Eq. (^), becomes 

The equation for the horizon area can be expressed as 

^=e""A^. (52) 

dv 

Given an explicit form for T^,„(t>), we need to integrate this pair of equations subject 
to the final conditions that 6* = and A = Af in the future. 

The 5-function pulse example presented in the Appendix indicates that the horizon 
tends to anticipate the energy flux by a time scale of order k~^ = AM. However, as 
shown there, the exponential factor e'^'" tends to strongly suppress any influence on 
timescales long compared to k~^. If we send in a negative energy pulse with duration 
much greater than M, followed by a positive energy pulse, we would expect that A 
will first decrease, and then increase, always a time ~ n~^ ahead of the energy pulse. 

4.2 Compact Pulses 

In this section, we construct a specific waveform for on the horizon which results 
in a (quasi-periodic) negative energy flux down the black hole. In general, for wave 
propagation in a black hole spacetime we can write 

= 0(t,r*,0,$) = :^^^^r,^, (53) 

r 

where the YJm are the usual spherical harmonics. 

Choose the ingoing modes to be / = only. To have an energy pulse which is 
bounded in time, we will use the following form for on the horizon: 

v\v-T,f gTJt.r*) 

<P[fJ) = T^TJ 7^, 0<V<To, (54) 

^ ^ Tn^ 2MVa^ - - u, k j 
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and 0(f) = 0ift'<0ort'> Tq. We will take qt to be a sum of two sinusoidal 
functions, as this form is an example which produces a negative energy flux down the 
hole. In particular, let 

gx = Oi cos uiv + ao cos uqv 

= ^ (e^"^!^ + e"^"^^^) + — (e'"^"^ + 6"*'^°^) . (55) 

Note that there would be a potential problem with simply choosing gT(t,r*) on the 
horizon to be the form given by Eq. (|55|) and "chopping" it; we could introduce 5- 
function terms or their derivatives into T^^. To avoid this, we chose the form of 
given by Eq. (^if), so that both and 0,^, go to zero smoothly at f = and v = Tq. 

We can choose 1/M ^ c<Ji > cuq, to get a form for gx that is approximately 
a constant + a sinusoidal part with a frequency much less than 1/M. Now let 
To = 271 /uq and Ti = 27i/uJi, where we choose Tq > Ti 3> M. The latter condition 
will insure that the time periods over which the flux is negative are long compared 
to M. This allows the changes in the Schwarzschild geometry to be quasistatic. 

The bounds on oq and ai are now restricted by the requirement that (1— Stt^ 0^) > 
everywhere. Our pulse is a wavepacket which consists of a range of both low and high 
frequency modes. The effect of the potential barrier is to reduce the value of on the 
horizon compared to its value just outside the potential barrier. The low frequency 
modes are the ones most attenuated by the potential barrier, whereas the high fre- 
quency modes are the ones least affected. Since increases roughly monotonically 
with decreasing r, the maximum value of would occur just outside the barrier, at 
r = To ~ 4M . On the other side of the barrier, the width of our pulse is ~ Tq; hence 
the wavepacket on the horizon will be dominated by modes with frequency ~ ujq. 
If we trace these modes backwards through the potential barrier, they will grow in 
amplitude by a factor of about iT^ooT ~ (4Ma;o)~ where T^^o is the transmission 
coefficient for a spherically symmetric wave of frequency uq- It is given by [^ 



%,o = -^iMuo, (56) 

for uq <^ 1/M. Note that we need both of the frequencies, ujq and Ui, to be nonzero 
so that the transmission coefficients for both components be nonzero. For this reason, 
the choice, Eq. (^, used in flat spacetime is not appropriate here. We will choose 
^ > 0, flo > oi, and ui ^ ujq, so 

■' '- "" (57) 



GAMuo V47r tq 

where we have used the fact that the factor of v'^{v — Tq)'^/Tq attains its maximum 
value of 1/16 at f = To/2. The requirement that (1 — 8tt^ (pmax'^) > gives 

|0max|<^==^, (58) 
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so therefore we have that 

2567rM2 [2 , , 

ao<^^^^, (59) 

where we have used cuq = '2^11 /Tq and tq ~ 4M. 

4.3 Two Illustrative Examples 

We now present two interesting numerical examples. The first is a case where A 
never decreases even though the flux is periodically negative. In this case, the GSL is 
rescued by the "precognition" of the horizon. That is, the horizon starts to respond 
acausally in anticipation of the subsequent energy pulse. This counter-intuitive be- 
havior of the horizon is discussed in more detail in the Appendix. In the second 
example, the periods of negative energy flux last long enough that this precognition 
effect cannot keep the horizon area from decreasing periodically, although the inte- 
grated area change is positive. What keeps the GSL from being violated during the 
periods of area decrease is the effect of the cfP' term in the formula for the entropy. 
These examples demonstrate that consideration of both effects is important for main- 
taining the GSL. If one or the other of them is left out, one can find examples where 
the law is violated. 

Let us choose units in which M = 1, so that k = 1/4. Strictly, this is the final 
mass, but here |AM|/M <^ 1, so we can assume that n is approximately constant. 
Then our final conditions in v are that at f = Tq, ^ = and A = Af = IQn. 

Case (i): Choose Tq = 40, Ti = 4, and ^ = 1/6. From Eq. (0), the bound on qq 



in our new units is Oq < (2567r/To) y 2/^, so in this case ao < 69.6. We conservatively 
choose ao = 0.16 and ai = 0.5ao. The behavior of the flux, the horizon area, and 
the entropy are represented in Fig. 1. The three quantities have been rescaled in 
order to display all three on the same plot. The plotted rescaled quantities are: 

Fscaled = 3.2 X lO^T.^, Ascaled = 10^ {A - Af)/Af, and Sscaled = lO^S - Sf)/Sf, 

where Sf = Af/A = An is the final entropy in our units. In this case there are three 
negative energy pulses embedded in the overall flux, but the duration of each pulse 
is short enough that the negative and positive parts of the flux occur close together 
in time. The acausal expansion of the horizon, in anticipation of the (larger) positive 
parts of the flux, is enough to save the GSL in this case. It is very important to 
take this effect into account in the analysis of the dynamical case. For a static black 
hole, the area of the horizon is simply IGttAP. This is also the area of the apparent 
horizon in the dynamical case. If one naively adopts this form for the area during the 
dynamical phase (as we mistakenly did in an earlier analysis of this problem ||31||), 
then even taking into account the change in the mass, one can still find cases where 
the GSL appears to be violated. 

Case (ii): Choose Tq = 400, Ti = 40, and ^ = 1/6. From Eq. (|5|), Oq < 6.96. 
We again choose Oq = 0.16 and ai = O.Soq. The behavior of the flux, the horizon 
area, and the entropy are shown in Fig. 2. The plotted rescaled quantities are: 

Fscaled = 1.6 X lO^T™, Ascaled = lO^A - Af)/Af, and Sscaled = lO' {S - Sf)/Sf. 
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Figure 1: Flux, horizon area, and black hole entropy as functions of v for Tq = 40, 
Ti = 4, ao = 0.16 = 2ai, and C, = 1/6. The plotted rescaled quantities are: Fscaied = 
3.2 X lO^T,,, Ascaied = 10^ (A - Af)/Af, and Sscaied = 10^ (S - Sf)/Sf, respectively. 
The behavior of the horizon area and the entropy are almost the same. 
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Figure 2: Flux, horizon area, and black hole entropy as functions of v for Tq = 400, 
T\ = 40, ao = 0.16 = 2ai, and C, = 1/6. The plotted rescaled quantities are: 

Fscaled = 1.6 X lO^T^^, Ascaled = 10^ (A - Af)/Af, and S.caled = lO^S ~ Sf)/Sf, 

respectively. Here the horizon area undergoes periodic decreases, but the entropy is 
always non-decreasing. 

There are again three negative energy pulses. However, in this case, each pulse lasts 
sufficiently long enough that the acausal behavior of the horizon is not sufficient to 
rescue the GSL, as evidenced by the three periodic decreases of the horizon area. 
That is, if the entropy were simply given by S* = A/4, as in the usual case, then the 
GSL would be violated during the periods of area decrease. What prevents this from 
occurring in this case is the influence of the scalar field term in the formula for the 
black hole entropy. Note also the phase lag between the change of sign of the flux and 
the onset of area decrease or increase. The scaling factors used in the two figures are 
different. The values of the scalar field on the horizon are approximately the same 
in both cases, and hence the magnitude of the fractional difference between S and A 
are of the same order of magnitude. However, the magnitude of the flux is smaller in 
Fig. 2 than in Fig. 1, because of the slower rate of change of 0. 
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5 Conclusions 

We have demonstrated that classical inassless non-minimally coupled scalar fields can 
be used to produce disturbingly large negative energy fiuxes, even in fiat spacetime. 
These classical negative energy fiuxes are temporary in the sense that the time- 
integrated fiux is always positive. However, unlike the situation in quantum theory, 
there are no constraints here analogous to the quantum inequalities. Therefore, it 
is possible for such fiuxes to be sustained over macroscopically long time intervals. 
These fiuxes appear to have magnitudes sufficient to produce gross violations of the 
second law of thermodynamics. 

However, whether such violations occur could depend on the details of how these 
fields interact with matter. In an attempt to circumvent this issue, we examined these 
negative energy fiuxes in the presence of a relatively unambiguous energy detector, 
a black hole. We considered scenarios in which a negative energy fiux, composed 
of classical non-minimally coupled scalar fields, is injected into the black hole. It 
was proved that the GSL is in fact always satisfied. This conclusion is true for any 
positive or negative value of the coupling parameter, subject to the condition that 
(1 — 87r,^0^) > 0, and is independent of the scalar potential V{4>). The proof is also 
not restricted to Schwarzschild black holes, and therefore will hold for rotating black 
holes as well. The validity of the GSL hinges on two key effects - the acausal behavior 
of the event horizon and the presence of an additional term in the formula for the 
black hole entropy which depends on the scalar field. 

In the dynamical case, the area of the event horizon, in contrast to that of the 
apparent horizon, depends upon the future history of the spacetime. This behavior is 
reviewed and illustrated in the Appendix. For a Schwarzschild black hole, when the 
temporal separation of the negative and positive energy pulses is no more than several 
times M, as was the case in Fig. 1, this effect alone can be sufficient to preserve the 
GSL. Here the horizon moves outward in anticipation of the (larger) positive energy 
in such a way that A is non-decreasing. 

If, however, the temporal separation becomes long compared to M, this effect 
is no longer sufficient to prevent A from decreasing. This behavior is illustrated in 



Fig. 2. A similar result was found in the case of Brans-Dicke theory in Ref. [^ 
Here the 0-dependence of 5* is crucial for insuring that the entropy is non-decreasing. 
Note that this 0-dependence is also not always sufficient by itself to save the GSL 
without the acausal behavior of the horizon. For example, the product of the area of 
the apparent horizon times (1 — Stt^^^) can decrease. 

There are still some questions left unanswered at this point. One is whether 
negative energy could actually destroy a black hole. Our treatment in Sect. |3.2| and 



those in Refs. p6| , p^ have assumed either a final black hole or cosmic censorship, 
and hence do not rule out this possibility. Another question is whether one could 
violate the GSL if (1 — Svr^^^) is not positive everywhere; this situation would involve 
either super-Planckian values of the scalar field or very large values of the coupling 
parameter. It might also result in a sign change of the effective local Newton's 



constant, given by Eq. (|3T|), and non-positivity of the entropy density. Recall that 
this condition played a crucial role at several steps. We assumed it in our proof of 
the GSL, and it is also needed to define the transformations from the Jordan to the 
Einstein frame, and hence to make the argument that the area in the Einstein frame 
provides a suitable non-decreasing quantity. It may also be necessary for the initial 
value problem to be well-posed [^. If (1 — Svr^^^) — > 0, then one would expect 
the effective stress tensor to become singular. This should lead to a backreaction 
effect which will prevent the system from actually reaching (1 — Stt^*^^) = 0; such an 
effect was indeed found in the case of cosmological spacetimes in Ref. |Q . However, 
T^u can remain finite if the numerator in Eq. (^ vanishes at the same points that 
(1 — 8TT^(f)'^) —>■ 0. This is what happens in some of the solutions of Barcelo and 
Visser |l^, |16|, [l^. Whether one could similarly construct examples which violate 
the GSL when (1 — 87r^0^) < remains unknown, although we strongly suspect so. 
Nonetheless, it is quite interesting that, even for relatively weak fields and coupling 
constants: (a) one can get large sustainable, albeit temporary, classical negative 
energy fluxes, and (b) that such fluxes do not violate the GSL. 

The acausal nature of the area of the event horizon is a disturbing feature of 
the current formulation of black hole thermodynamics, and one which has drawn the 



attention of various authors [35|. It is not clear whether this feature will survive in 
future theories. A definition of entropy based on a statistical mechanical enumeration 
of states would seem to have to be causal. This is a deep question which may have 
to await a more complete quantum theory of gravity for a resolution. 

Finally, there is the question of whether the large negative energy associated 
with a classical scalar field could produce any dramatic effects, such as violations of 
the second law for an ordinary thermodynamic system. More generally, does this 
negative energy produce any observable effect on ordinary matter? It would seem 
rather peculiar if one could violate the ordinary second law but not the GSL. One 
would suspect that they should either both hold or both be violated. However at 
present we do not know of a similar general argument that would guarantee that the 
ordinary second law is not violated. 

If such classical scalar fields can exist, that raises the question of why only cer- 
tain classical fields are capable of producing large energy condition violations. For 
example, it is often thought that the conformally coupled scalar field is "physically 
reasonable" because it faithfully mimics certain aspects of the electromagnetic field. 
But why then does the classical electromagnetic field obey the WEC, while the con- 
formally coupled scalar field does not? 

It is usually easier to get negative energy in the context of quantum rather than 
classical fields. Moreover, in the regime of quantum fields, there seem to be some 
rather strong restrictions imposed on the magnitude and extent of negative energy 
densities and fluxes, in the form of the quantum inequalities. Given that nature seems 
to tightly restrict quantum violations of the weak energy condition, how seriously 
should one take classical violations? These are all issues for further enquiry. 
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6 Appendix 

In this appendix we analyze the simple case of a spherical (5-function pulse of positive 
energy imploding onto an already-existing black hole. We examine this case for two 
reasons. The first is to illustrate the somewhat counter-intuitive dynamical behavior 
of the horizon which is important in our discussion. The second is to present this 
calculation in a pedagogically detailed form - one which we have been unable to find 
in the standard literature. 

For simplicity, we consider a Schwarzschild black hole of initial mass Mq, which 
absorbs a spherically symmetric (5-function pulse of positive energy null fluid and 
becomes a black hole of mass M. Let A be the afiine parameter on the horizon and 
suppose that the pulse intersects the horizon at A = 0. In the case of a spherically 
symmetric pulse the Raychaudhuri equation for the null generators of the horizon 
reduces to the simple form 

^ = 4-V*^r. (60) 

For positive energy R^yk^k^ = SnT^^k^k^ > 0, so set R^yk^^k^ = a6{\), where a is a 
positive constant. 

Let us first solve the Raychaudhuri equation without the 5-function term: 

de e^ 

This equation has the simple solution 

After the absorption of the pulse the generators of the horizon have zero expansion, 
so we want ^ = for A > 0, so let 

''(^> = {2/(A + A„). A<0- «'3) 
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i. e. , ^ jumps discontinuously from 2/Ao to as A passes through 0. This gives us a 
(5-function term in d9/d\ of —25{X)/\q. Thus the above form of 6' is a solution of 



de 

dX 



o2 

-- - a6{X) 



(64) 



if a = 2/Ao, or Aq = 2/a. (The divergence of ^ at A = — 2/a represents a conjugate 
point in the distant past, which will be discussed later.) 

Let us now generalize our solution slightly by assuming that the pulse intersects 
the horizon not at A = 0, but at A = A+. Then our solution for 6 becomes 



^(A) 



, A > A+ 

2/(A - A+ + 2/a) , A < A+ 



(65) 



For the (5-function pulse case, in the region where A < A+, we can then integrate 
Eq. (!§) to get 

In A = 2 In (A - A+ + 2/a) + C . 



We fix the constant of integration C by the requirement that A 
area as A — > A+, so we have 

C = In A/ -2 In (2/a). 

Our expression for A then becomes 

'A-A+ + 2/ay 



or 



A{X) = Aj 



A{v) = Aj 



2/a 



, A < A+ , 



(66) 
Af = final horizon 

(67) 



(68) 



-jKV-l- 



+ 2k) 



2k 



f < f + 



where we have used Eq. (^Sj). We can write this in the form 



A{v) 



A, 



A/ { [a(e'^(''-^+) -l) + 2/te-™+] /{2k 6''^''+)] 



V >v+ 



(69) 



(70) 



If |f — f+l > K ^^ and since f — f+ < 0, then e"^" ^+^ ~ 0, and we have 

,2 



A 



A 



f 



2k 6-^"+ -a 



2k €-""+ 
a 



AAl-—e 
2k 



KV+ 



Ao. 



(71) 



Thus on a timescale of order k~^ = AM, A increases from Aq to Af. In the distant 
past (corresponding to A — >• A+ — 2/a), A ^ 0, as can be seen from Eq. (pSD, but 
this is the conjugate point associated with the formation of the black hole or with 
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the past singularity in the case of an eternal black hole. After the black hole forms, 
its area remains nearly constant until a few times M before the positive pulse arrives 
at f = w+. 

We see that the horizon "anticipates" the subsequent arrival of the pulse. This 
counterintuitive acausal behavior arises from the way the horizon is defined. Recall 
that the horizon is defined as the boundary of past null infinity, i. e., the boundary 
of the region from which it is possible for light rays to escape to infinity (see, for 
example, Ref. PB[, p. 300). Thus the horizon has the peculiar property that one 



must know the entire future evolution of the spacetime in order to know where the 
horizon intersects any given spacelike slice. For a static black hole, the null generators 
of the horizon have zero expansion. In order for the null generators in the horizon of 
the final (static) black hole to have zero expansion after the absorption of the pulse, 
the generators must have nonzero expansion prior to the pulse's arrival. The null 
rays in the horizon of the final black hole are rays that would have escaped to infinity 
were it not for the arrival of the positive pulse, which refocused them to have zero 
expansion. (See Fig. 59 of Ref. [Q, and the discussion in Ref. ||37[| .) 

Our calculation assumes that k is approximately constant and hence applies when 
the increase in area is small: 

^ Ao(l + - 6^"+ 
\ K 

= Ao + AA. (72) 

If we take Aq = IGvrMo^ to be the initial area of the black hole in the distant past, 
where Mq is its initial mass, then we have 

AA a ,„ AM , , 

^ - e^"+ = 2 -— . 73 

Aq K Mq ^ ' 

In this approximation, the change in the mass of the black hole is 

AM = —Mo 6^"+ 
2k 



= ^e-, (74) 

where we have used k = 1/(4M) ^ l/(4Mo). This agrees with the result obtained 
by calculating the change in mass directly from Eq. (p2|) . 

Let us now check that the conjugate point which appears in our solution occurs 
in the distant past. From Eq. ( |65D we see that the conjugate point occurs at 

A = Ae = A+ - - . (75) 

a 
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Using Eqs. (^) and ([7^, one can obtain 

2 AM 

a = 

MoA+ 

and so 



(76) 



+ AM ^ V AM/ ^ ^ 

Since in our approximation AM <^ Mq, Ac < and |Ac| ^ A+. Note that A = 0, or 
equivalently v = — cxd or \^ = 0, corresponds to the intersection of the past and future 
horizons in the maximally extended Schwarzschild spacetime. Thus for an eternal 
black hole, the conjugate point cannot occur to the future of the past horizon. In 
the case of a black hole formed from collapse, the same parameterization of A yields 
A > everywhere on the future horizon outside of the collapsing body. In either case, 
Ac < implies that the conjugate point occurs in the distant past. 
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